Improved Light Absorption by Quantum Confinement and Band 
Folding: Enhanced Efficiency in Silicon Based Solar Cells 
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The improvement of light absorption in Si/BeSeo.4iTeo.59 heterostructures for solar cell applica- 
tions is studied theoretically. First, using simple approaches we found that light absorption could be 
improved in a single (uncoupled) quantum well with a thickness up to 20 A. Second, by semiempir- 
ical tight-binding methods we calculated the electronic structure and optical properties of various 
(Si2)n/(BeSeo.4iTeo.59)m [001] superlattices. Two bands of interface states were found in the band 
gap of bulk Si. Our calculations indicate that the optical edges are close to the fundamental band 
gap of bulk Si and the transitions are optically allowed. 

PACS numbers: 73.21. Cd,73.21.Fg,78.67.De 



I. INTRODUCTION 

Silicon is the principal material for integrated circuits. 
Not only is the band gap (1.12 eV) ideal for room tem- 
perature operation, but also the oxide (Si02) provides 
the necessary flexibility to fabricate millions of devices 
on a single chip. High integration implies high-speed op- 
eration that is limited by the interconnect propagation 
delay of the signal between devices. This constraint sug- 
gests that the integration of Si micro-electronics might be 
aided by optical interconnection. Unfortunately, silicon 
does not respond strongly to optical exitations because it 
is an indirect band gap semiconductor: the band extrema 
for electrons and holes are located at different points in 
the Brillouin zone (Fig. |l|). Therefore, intrinsic formation 
or recombination of electron-hole pairs becomes a three- 
particle event, which is weaker than a two-body process. 

There has been a tremendous effort in exploring ways 
of breaking the silicon lattice symmetry and mixing 
diffcijent momentum states in order to induce optical 
gairala. The radiative efficiency depends on the compe- 
tition between non-radiative fast processes and relatively 
slow radiative processes. To optimize the efficiency we 
have to eliminate non-radiative channels by having high 
purity materials and increase the oscillator strength of 
radiative channels. The above criteria also apply to pho- 
tovoltaics. 

One modality for improving the optical response is 
quantum confinement. Confinement of the charge density 
in quantum wells (QWs) permits the relaxation of the op- 
tical selection rules for interband transition. In addition, 
the band folding in superlattice (SL) structures will en- 
hance the absorption. Consequently a Si based SL will 
reduce the symmetry which translates into band folding 
toward the zone center; and as a result vertical transi- 
tions will be available at energies closer to the indirect 
band gap. 

Beryllium-chalcogenides are good candidates for Si- 
based heterostructures. They are wide-band gap zinc 
blende semiconductors with lattice constants close to 
that of Si. Thus BeTe and Be Se have the lattice con- 




FIG. 1: T 
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band structure of bulk Si with the parameters of 
adapted with the spin-orbit constant A = 0.045 
eV. The single group notations are used since the spin-orbit 
splitting cannot be distinguished on the chosen scale 



stants of 5.6269 and 5.1477 A, respectively, 3.6 % larger 
and 5.2 % smaller than Si. Vegard's law indicates that 
the lattice matched composition with Si is BeSeg. 41X69.59. 
Recent developmcntfflcl in the growth of silicon lattice- 
matched BeSeo.4iTeo.59 open the opportunity for a new 
class of Si based devices. 

We analyze quantum confinement and band folding 
in silicon based heterostructures (Si/BeSeTe heterostruc- 
tures) as ways of enhancing light absorption and, there- 
fore, improving the efficiency of solar cells. In this study 
we use a semiempirical tight-binding (TB) methodQ. In 
section II we present the tight binding parameters of 
BeTe and BeSe, fitted to fi^^t principles calculations. Thg 
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model is that of Vogl et al, 3 modified according to ChadiQ 
as to include the spin-orbit interaction. The bulk TB 
parameters of BeSeo.41Teo.59 are considered as virtual- 
crystal approximation of constituent materials. The TB 
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parameters of silicon are those of Vogl et al.p l adapted to 
include the spin-orbit coupling. In section III we discuss 
some methods to improve the light absorption in Si based 
heterostructures such as quantum confinement in a sin- 
gle (uncoupled) QW and band folding in SL structures. 
A significant increase in direct absorption is obtained for 
a narrow uncoupled silicon QW. However, thin barriers 
are also required in order for carriers to move across the 
structure. In this spirit, we calculate the electronic and 
optical properties of Si/BeSeTe SLs using the mentioned 
TB method. In the last section the conclusions are out- 
lined. 



II. NEAREST NEIGHBOR TIGHT-BINDING 
PARAMETERS OF BERYLLIUM 
CHALCOGENIDES 

Since both Si and BeSeo.41Teo.59 are indirect band ma- 
terials, full band calculations are required for the theo- 
retical understanding of electron transport and optical 
properties in such heterostructures. For this purpose we 
use the empirical tight-binding (ETB) method as one of 
the most used tool£,in research of complex molecular and 
solid-state systemaj. Despite the fact that ETB is based 
on physical approximations such as a one-particle pic- 
ture, short-range interactions, etc. (aposteriori justified, 
however), it gives fast and satisfactory results. In terms 
of accuracy, it lies between the very accurate, very expen- 
sive ab-initio and the fast but limited-accuracy empiri- 
cal methods. TB techniques are increasingly employed 
in structures like resonant tunneling diodes, QWs, and 
SLs. They are more complete than envelope-function or 
k-p methods because they incorporate the entire band 
structure of the constituent bulk materials in a transpar- 
ent manner. Also the short-range nature of the model is 
suitable for modeling heterointerfaces which are present 
in such quantum structures. In ETB models the elec- 
tronic wave functions are expanded as a linear combina- 
tion of atomic orbitals. The real attractiveness resides in 
the great versatility of the method that can be used for 
a variety of physical situations like: (i) molecular mod- 
els based on linear combination of atomic orbitals with 
strong directional properties; (ii) calculation of electronic 
states in unknown systems based on transferability of 
the Hamiltonian matrix elements as determined in known 
cases (i.e. transferability of the Hamiltonian matrix ele- 
ments from bulk semiconductors to low dimensional sys- 
tems like QWs, quantum dots, etc.); (iii) quantitative 
predictions of physical properties (like fitting the band 
structure of a semiconductor to ab-initio and experimen- 
tal data and using those parameters to calculate opti- 
cal properties of the same semiconductor). The last two 
points will be used in this paper. 

0ns, of the most popular TB models is the sp^s* 
modefi. This iS|-a twenty-band model if spin-orbit cou- 
pling is includedQ. The presence of additional s* orbital 
is able to reproduce the band gap in indirect semiconduc- 



TABLE I: Matrix elements in eV of nearest neighbor sp^s* 
model incl uding spin-a rbit interaction. The notation is ac- 
cording to Vogl et aLp . Virtual-crystal approximation was 
used to calculate the matrix elements of BeSeo.41Teo.59. 
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BeTe and BeSe are quite new materials in the sense 
that there are few experimental facts about these semi- 
conductors. The bulk TB parameters of BeTe and 
BeSe are determined by fitting the ab-initio calcula- 
tions of Fleszar and Hankecl. Virtual-crystal approxima- 
tion (VGA) is used to determine the TB parameters of 
BeSeo.4iTeo.59 . The approximation is expected to work 
well because BeTe and BeSe crystallize in zinc-blende 
structures and have similar band structures. The nearest 
neighbor parametrization in the sp^s* model requires 15 
parameters including spin-orbit interaction. The spin- 
orbit coupling was included due to the large spin-orbit 
splitting of the valence band, of about 0.45 eV for BeSe 
and 0.96 eV for BeTe. To find the nearest neighbor 
parametrization we used a similar procedure to Ref. 
The procedure is more consistent than that presented in 
Ref. H We adapted the procedure in order to consider 
spin-orbit coupling. The results are shown in Table |. 

The focus of these calculations was to find the best 
parameters that reproduce the valence band edges (Fg 
and F7 and conduction band edge at Fg and Xi). The 
reproduction of band edges at the zone center of the va- 
lence and conduction bands and the band gap between 
valence band and conduction band is within 3% error. 
The spin-orbit coupling for Be was chosen to vanish, be- 
cause Be is a light element and its spin-orbit coupling 
is negligible. This is also consistent with the spin-orbit 
splitting of 0.1 eV for another Bc-chalcogenide, BeS with 
the sulfur atom, lighter than selenium and tellurium 
atoms. The energy band diagram for BeSeo.41Teo.59 is 
shown in Fig. |^. The material appears to be indirect 
band gap with the conduction band minimum at the 
X point. The predicted band gap is 2.97 eV. Previ- 
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of the Si/BeSeo.4iTeo.59 SL. 




FIG. 2: Energy bands for BeSeo.41Teo.59 calculated with TB 
parameters from Table 1. 

ous measurementsi indicate a conduction band offset of 
1.2 eV for the Si/BeSeo. 41X60.59 heterostructure. This 
leaves a valence band offset of about 0.65 eV for the 
Si/BeSeo.4iTeo.59 heterostructure. The TB parameters 
and the valence band offset will be used in the next sec- 
tion in order to calculate electronic and optical properties 



III. IMPROVED OPTICAL ABSORPTION BY 
QUANTUM CONFINEMENT AND BAND 
FOLDING 



A. A Single Quantum Well 

The absorption edges in bulk silicon are indirect and 
the jumping of the electron from valence band (VB) 
to conduction band (CB) is assisted by a phonon (in 
electron transitions from VB to CB the k momentum 
must be conserved). The band structure of Si is shown 
schematically in the Fig. |l|. The indirect band gap is 
Ai ~^25^' absorption process involves two steps, one 
is electron-photon interaction and the other is electron- 
phonon interaction. Basically we have two processes de- 
noted by 1 and 2. In process 1 an electron is first excited 
to the Fis state and then by an emission of a phonon the 
electron arrives in the Ai CB state. Similarly for process 
2, however, in this process it starts with an emission of 
a phonon. Mathematically, the transition probability is 
given by Fermi's Golden Rule for second order perturba- 
tion 
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Eio - hw 



S [Ec (kc) - Ev (kv) ~ hw ± htOp) , 



(1) 



where jO) represents the initial state, \i) the intermediate 
state and|/) the final state. Hep and HeR are the Hamil- 
tonians for electron-phonon and electron-photon interac- 
tion, respectively. The energies hw and fiiOpave the ener- 
gies for photons and phonons, respectively. The second 
sum on the right-side hand is the oscillator strength of 
the transition. For this reason, there is a long absorption 
tail between 1.12 eV and and about 3 eV that refiects the 
indirect nature of the band gap. The sharp rise in ab- 
sorption with increasing photon energy starting around 



3.2 eV (380 nm) is associated with the direct transition 



at F point (F 



25' 



F15) whose energy is 3.4 eV (365 



In direct band gap materials light absorption in 2D 
systems is formally similar to 3D systems. Although the 
light absorption in indirect band gap 2D systems looks 
to be similar to 3D systems, this is not the case. The 
absorption coefficient for an uncouplod. Si/ BeSeo.41Teo.59 
QW has basically three componcntsllil 



J 



pij {Eg + Ee+ Eh - Tiw) + nx, S {Eg + E^ + Eh - huj - hw^ 



eh.Xa 



E P^eh ("Ag + 1) S {Eg +Ee+ Eh - huj + hwxq) 



eh.Xn 



(2) 



where ^4 is a constant, gives the direct bandgap con- with 1-phonon absorption, andp^^ is the phonon assisted 
tribution and p"^ gives the phonon assisted contribution 
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FIG. 3: The overlap integral (Eq.(|5f) for infinite wells in sil- 
icon. The wave vector ko in Eq.(p|) is the wave vector for 
silicon conduction valleys along z— axis. 



contribution with 1-phonon emission. Basically, 



Peh = Pcvleh (ko) , (3) 



PthM = Pth,\q = Pcvleh (q - ko) i?A, (4) 

Pcv is the bulk dipole matrix element between bands, Rx 
is the matrix element contribution from electron-phonon 
interaction, and leh is the overlap between the envelope- 
functions at different Brillouin points: 

leh (q) = Jdr V: (r) (r) 6"*°^ (5) 

where ko is the location in fc— space of each conduction 
valley. For silicon, if z is the growth direction, then the 2 
valleys along z-directions (with ky = = 0) are respon- 
sible for direct transitions. The other four valleys con- 
tribute to phonon-assisted absorption. For infinite wells 
the overlap integral /e^ has the behavior depicted in Fig. 
|. For very narrow QW's the overlap integral tends to 
one, which means that in a genuine 2D system only the 
transverse momentum has to be conserved. For the other 
asymptotic limit, i.e. very wide well, the overlap integral 
vanishes. In physical terms, this says that the system 
became genuinely 3D and any electromagnetic transition 
has to be vertical in the absence of phonons. Moreover, 
from Fig. ^ we may expect to have a strong direct tran- 
sition (with an overlap integral no less than 0.1) for a 
QW with a width up to 20 A, i.e. an ultrathin QW. In 
the following we will discuss the possibility of improving 
the oscillator strength in Si/ BeSeo. 41X69.59 based SLs. 



B. Electronic Structures of (Si2)n/(BeSeo.4iTeo.59)m 
[001] Superlattices 

The use of SL structures relaxes the condition of ul- 
trathin QW. However, the SL has to have thin barriers. 
Thin barriers allow photo-generated carriers to be swept 
away by the built in electric field without recombination. 
Thicker barriers induce localization of carriers and im- 
plicitly higher recombination rate. The alternation of 
QWs and barriers along [001] growth axis will generate 
a band folding in the^SIi-and band mixing of zone-center 
and zone edge statesE3Ea. The band mixing will inher- 
ently enhance the oscillator strength for direct transition 
in Si structures. On the other hand, band folding will 
induce states to which vertical transitions are possible at 
energies by far lower than 3.4 eV, the lowest energy for 
vertical transitions in bulk Si. 

We consider a Si/BeSeo.4iTeo.59 SL whose layers 
are perpendicular to [001] direction. We employ 
the nearest neighbor sp^s* TB Hamiltonian includ- 
ing spin-orbit interaction. We denote this SL as a 
(Si2)n/(BeSeo.4iTeo.59)m SL with n two-atom thick lay- 
ers of Si and m two-atom thick layers of BeSeo.41Teo.59 
repeated periodically. The defined supercell consists of 

2 (n + m) adjacently bonded atoms Si, Si, Si, Be, 

Te/Se, Be, Te/Se, Be, Te/Se. The TB matrix el- 
ements of the Si/BeSeo.4iTeo.59 SL are taken over di- 
rectly from those bulk values. The on-site energies of 
the BeSeo.4iTeo.59 are accordingly changed to match 
the valence band offset at the interface. Simple av- 
erages were used to supply the parameters connecting 
different materials at the interface. Since the spin or- 
bit must be included the SL Hamiltonian will be rep- 
resented as having 20 {n + m) functions. Once the TB 
matrix elements have been established the SL band struc- 
ture reduces to the diagonalization of the 20 {n + m) x 
20 {n + to) matrix Hamiltonian. The band structure of 
(Si2)n/(BeSeo.4iTeo.59)m SL for m = 71 = 2, m = rt = 4, 
m — n — 6, and m = n = 8 are displayed in Fig. ^. The 
zero of energy corresponds to the top of the valence band 
in bulk Si. The band folding effect can be seen as many 
crowded subbands. Two interface bands {ly and Ic), one 
empty and one occupied, were found. They lie in the 
lower and upper parts of the band gap of bulk silicon, 
respectively. The origin of these interface bands rests 
on the polar nature of the interface as was also found 
in GaAs/Ge SLIIj. The polarity of the interface origi- 
nates from the large differences in the on-site energies 
for the constituent atoms (Si and Be or Se/Te). Even if 
a (110) non-polar interface is used, one interface band is 
still found in H-VI/IV SLsBEa. We calculated the planar 
charge density of some of the occupied and empty band 
edge states of the SL with m = n = 8 at F point in the 
Brillouin zone. The planar charge density of some of the 
band edge states and interface states are depicted in Fig. 
^ for the SL with m = n = 8. By interface states we 
mean states which die away within few layers from the 
interface. We denote by F^", the interface states of the 
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ly band and by F^^, the interface states of the Ic band 
at r point. The charge density of the T^"" interface state 
has a maximum at the Be-Si interface, while the next 
occupied state is confined, but stiU is locahzed toward 
the Si-Se/Te interface. Only the third occupied state is 
genuinely confined in silicon slab. The charge density of 
the r^'^ interface state has a maximum at the Si-Se/Te 
interface, while the next empty state is confined, but still 
is localized toward the Si-Be interface. Again, only the 
third empty state is genuinely confined in silicon slab. 
We believe that this is associated with acceptor behavior 
oij-Be and donor behavior of Se and Te with respect to 



C. Optical Properties of (Si2)n/(BeSeo.4iTeo.59)m 
[001] Superlattices 

The electronic contribution to absorption spectrum is 
given by aabs (w) '^rWe2 (w), whereo-Qbs {^^) is the absorp- 
tion coefficient and@ 

£2 (w) = fcv,k5 {Ec.k - Ey^k - fuv) (6) 

muj\l ^ — ' 

c.v.k 

is the imaginary part of dielectric function. Here m is the 
electron mass, D, is the volume, e is the electron charge, h 
is the Planck constant and /^.^is the oscillator strength 
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FIG. 5: Planar charge densities in the (Si2)8/(BeSeo.4iTeo.59)8 [001] superlattice at F point, (a) Occupied states: the interface 
state F^*^, and the first top two confined states (F^^ and F^^). (b) Empty states: the interface state F^"", and the first two 
confined states (F'^^ and F'^^). The solid vertical lines denote the interfaces. 



for the direct transition from the state |w,k) to |c, k), 
with the photon momentum neglected. The oscillator 
strength is defined as 



energies. Similar results were found in the calculations 
for porous Si with periodic boundary conditiontcil. Eq. 



/c 



2 |(c,k| £-p \v,k) 



Er 



v.k 



(7) 



In Eq. (|^) |w,k) and |c, k) are the valence and con- 
duction band, eigenstates, E^^k and Ec,k are their corre- 
sponding energies, e is the polarization of light, and p is 
the momentum operator. In the empirical tight-bindips 
approach the momentum matrix element is defined asllfl 



(c,k|p|c,k) = -(c,k| Vki?(k) |i;,k) 



(8) 



Because the variation oLthe oscillator strength over 
the Brillouin zone is smal£3, we first calculate the joint 
densities of states (JDOS). We assume that light is prop- 
agating along the SL growth direction. The JDOS repre- 
sents the number of states that can undergo energy and 
k-conserving transitions for photon frequencies between 
uj and uo + duj. The JDOS associated with Eq. (^) are 
shown in Fig. |6| for SLs with m = n = 2, m = n = 4, 
m = n = 6, and m = n = 8. A 0.05 eV broadening was 
considered for each electronic energy. The summation 
over Brillouin zone was replaced by-the summation over 
special points in the Brillouin zonellSE^I. 

Due to band folding the absorption edges for vertical 
transitions are lowered toward the indirect band gap of 
bulk Si. Moreover, the curves rise slowly and, with in- 
creasing m (= n), the absorption edges extend to lower 




Energy (eV) 



FIG. 6: Joint density of states corresponding to optical tran- 
sitions for (Si2)n/(BeSeo.4iTeo.59)m [001] superlattices with 
n = m = 2,4, 6,8. 
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FIG. 7: The oscillator strengths / of several interband tran- 
sitions relative to the direct transition in bulk Si (rjg' 
Tis) fbuikSi for the (Si2)4/(BeSeo.4iTeo.59)4 [001] superlat- 
tice. 

tells us that the strength of the optical absorption 
is also determined by the oscillator strength. We denote 
by VI, V2, and V3 the first three top valence subbands, 
and by CI, C2, and C3 the first three conduction sub- 
bands. We calculated the oscillator strengths of several 
interband transitions relative to the oscillator strength of 
the direct transition in bulk Si. The results are shown 
in Fig. for TO = n = 4. The oscillator strengths of 
interband transitions are at least 10 times smaller than 
their bulk counterpart and basically range from 10 
to 10 relative to T^^,^ — > Fis transition in bulk 
Si. The strongest transitions are those coming from the 



third confined hole level to the second and third confined 
electron levels. 



IV. CONCLUSIONS 

We studied light absorption of silicon-based het- 
erostructures (Si/BeSeo.4iTeo.59) that might be used 
in solar cell applications. We found that significant 
absorption occurs for thin uncoupled quantum wells 
with their width up to 20 A. We also studied the 
(Si)„/(BeSeo.4iTeo.59)m superlattice structures with a 
semiempirical tight-binding method. In order to per- 
form the superlattice calculations, we determined the 
tight-binding parameters for BeSe and BeTe in the sp^s* 
nearest neighbor model including spin-orbit interaction. 
Electronic structure and optical properties were calcu- 
lated for various superlattice structures. Two interface 
bands were found to exist in the band gap of bulk silicon. 
We found that band folding induces vertical transitions 
near the indirect band gap of bulk Si. In addition, cal- 
culated oscillator strengths for vertical transitions near 
the optical band edge show the mixing of the zone-center 
and zone edge states of bulk Si for conduction subbands. 
Therefore, the transitions are optically allowed and the 
response of silicon based heterostructures to illumination 
is enhanced. 
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